Photon statistics from a resonantly driven quantum dot 
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Photon correlations in the emission of a resonantly driven quantum dot are investigated, account- 
ing for the influence of the solid-state phonon environment. An analytical expression is derived for 
the second-order fluorescence intensity correlation function, from which regimes of correlated and 
uncorrelated photon emission are predicted as the driving field is varied. Experiments to investigate 
this effect would provide valuable insight into quantum dot carrier-phonon dynamics and are feasible 
with current technology. 
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There has been a great deal of progress over the past 
few years in refining fabrication and characterisation 
techniques for samples of small semiconductor quantum 
dots (QDs)^ This has lead to the ability to control and 
manipulate, with ever increasing precision, the quantum 
states of charge carriers confined within single QDs^ and 
QD arrays.— Such systems are prototype solid-state quan- 
tum processors, ideal for performing studies of coherent 
control techniques, decoherence, and quantum measure- 
ments that will be key in assessing the prospects for larger 
scale solid-state quantum information processing (QIP). 
QDs are also promising candidates for efficient single pho- 
ton sources,^ opening up potential applications in other 
paradigms such as linear optical- and hybrid light-matter 
QIP^ 

There is therefore considerable interest in employing 
QDs both as qubits themselves and as elements in a larger 
optical setup. For either application an understanding 
of the coherence properties of their optical transitions 
is key. In the former case to ensure high fidelity qubit 
operations, in the latter case to evaluate the degree of 
photon indistinguishability in the QD emission. Hence, 
an experiment that establishes both the nonclassical na- 
ture of the QD fluorescence, while at the same time ex- 
ploring the potential for coherent control of the QD ex- 
citonic states themselves, is desirable. In fact, a first 
experiment along exactly these lines has recently been 
reported by Muller et al.^ who measured the second- 
order fluorescence intensity correlation function, g^ 2 '(r), 
for a coherently-driven InGaAs QD, observing damped 
oscillations in the signal. 

Here, I shall show that not only does such a measure- 
ment provide a clear signature of a solid-state two-level 
quantum emitter, but that it is also a sensitive probe 
of the coherence properties of the QD excitonic ground- 
state. In contrast to the excitonic occupation dynam- 
ics usually investigated ) 9 ' 10 ! 11 ' 12 in the following I shall 
derive an analytical expression for the fluorescence in- 
tensity correlation function of a driven QD. I shall use 
it to show how the interplay between laser-induced co- 
herent exciton evolution, spontaneous photon emission, 
and dephasing due to exciton-phonon interactions 1 ^ man- 



ifests itself in the emitted photon statistics. In particular, 
it is expected that the phonon-induced dephasing rate 
should display a non-monotonic driving strength depen- 
dence giiHilliU This can be understood from a general 
resonance argument jSiiS the induced damping is dom- 
inated by phonons that are most strongly coupled to 
the QD carriers, which corresponds to those that have 
wavelengths comparable to the QD size d (or frequencies 
~ u/d, where u is the sound velocity^ -). These modes are 
most effectively driven when the carrier dynamics is reso- 
nant with them, which in turn implies that the Rabi fre- 
quency (QD-driving coupling strength) is resonant with 
them too. Hence, as the Rabi frequency is increased 
towards a value of u/d the carrier-phonon coupling in- 
creases as well, leading to a larger dephasing rate, only 
for it to decrease again as the Rabi frequency is increased 
further, moving the carrier dynamics out of resonance 
with the dominant phonon modes. I shall demonstrate 
below that these dephasing-rate variations give rise to a 
distinctive feature in the second-order correlation func- 
tion, namely the appearance, disappearance, and subse- 
quent reappearance of long-lived oscillations in g^(r) as 
the driving field is changed. 

We consider a single QD modelled by a two-level sys- 
tem with ground-state |0) (the semiconductor ground- 
state) and excited state \X) (a single exciton), separated 
by an energy ujq {% = 1). The dot is addressed by an 
external classical laser field of frequency ui and dot-laser 
coupling given in the dipole approximation by J7, with 
<C ujq. Both the phonon environment and the ra- 
diation field are modelled by harmonic oscillator baths 
of frequencies w q and #k, respectively, with q and k 
the corresponding wavevectors. The carrier-phonon in- 
teraction acts only on the dot excited state \X) and is 
taken to be of the usual spin-boson forma l 10 ' 11 ' 14 ' 15 ' 16 
Within a frame rotating at frequency loi this leads to the 
Hamiltonian>±£ H = v\X) (X\ + §(|0) (X\ + \X) (0|) + 
X) (X\ £_( ffq &t + 5q 6 q ) + £ k 6 k ata k + 



/k<7+ak), with rotating- wave approxima- 



tions on both the driving and system-photon coupling 
terms. Here, v = luq — u>i is the detuning of the laser from 
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the dot transition, <r_ = |0) (X\, a+ = \X) (0|, 5 q (/ k ) 
defines the exciton-phonon (exciton-photon) coupling, 6 q 
(ak) are the phonon (photon) annihilation operators, and 
an irrelevant term has been dropped. 

On taking a trace over the phonon modes, which are 
assumed to be in thermal equilibrium, the resulting spec- 
tral density J(ui) = J2q l5q| 2 ^( CJ — w q) completely de- 
scribes the interaction between the QD and phonons. 14 
For arsenide QDs, the coupling of the confined exci- 
ton to acoustic phonons by means of the deformation 
potential tends to dominate the dephasing dynamics, 
over the piezoelectric interaction or coupling to opti- 
cal phonons . 13 ! 15 ' 16 In this case, the coupling constant 
is given by g q = qD(q)/ v /2Aoj q F,— ^ where A is the 
sample density and V the unit cell volume. The form 
factor is £>(q) = / dr (D ft |Vft(r)| 2 - D e \^)\ 2 ) e" iq ' r , 
with tph and ip e being the confined hole and elec- 
tron ground-states, respectively, while Dh and D e 
are the corresponding bulk deformation potential con- 
stants. For clarity, I take a spherically symmet- 
ric harmonic confinement potential for the QD, giv- 
ing ip e {h) = (de(h)V^)~ 3/2 exp (-r 2 /2d 2 e{h} ), where d e{h ) 
is the electron (hole) ground-state localization length. 
For d e = dh = d we obtain J(u>) — u 3 (Dh — 

J 1 d 2 9 

D e ) 2 e 2^-/(47r 2 Aii 5 ) = aw 3 e~^/" c ) , where a linear 
dispersion w q = uq has been assumed. The spectral 
density is therefore super-ohmic with a natural high fre- 
quency cut-off at to r-j cj c — \[2ujd due to the finite QD 
size. 

Under resonant continuous-wave excitation the system 
reaches a quasi steady-state (ss)j^ after which time the 
normalised second-order fluorescence intensity correla- 
tion function of the QD, g^(r) = (I(t)I(t + rj) / (I(t)} 2 , 
may be written^ 

gW(r) = (a+a-)-' (l + (a z (r)) p(0 )=|o>(o|) • (1) 

This is proportional to the probability of detecting two 
emitted photons separated by a time r, and normalised 
by the form for independent detection. Here, and in the 
following, Oi (for i = x,y, z) are the usual Pauli matri- 
ces in the basis {|0), \X)}. Note that the evolution of 
a z (t) must be obtained conditional on the system being 
prepared in its ground-state at t = (p(0) = |0) (0|), 
i.e. detection of the first photon initializes the system 
state. The continuous-wave excitation condition ensures 
that there is no effect of any laser-pulse temporal profile 
in subsequent measurements. 

During their evolution the QD states are coupled both 
to a bath of phonons and the radiation field. Utilising 
the time-convolutionless projection operator technique 21 
a master equation for the reduced system density op- 
erator may be derived from the Hamiltonian and the 
Liouville-von Neumann equation, here within the Born 
approximation and up to second order in |</ q | and |/k|. 
Typically, the memory time within the radiation field is 
extremely short, r rn ~ (1/loq) ~ 1 fs for ujq of around 
1 eV. 22 Since fi <C ujq, we are then justified in treating 




12 3 4 

n (meV) 



FIG. 1: (Color online) D (red, solid line) and Re(^2) (yellow, 
dashed line) as a function of driving strength Q. Parameters: 
D e = -14.6 eV, D h = -4.8 eV, d = 4.5 nm, u = 5110 ms _1 , 
A = 5370 kgm" 3 , and T = 50 K, giving a = 0.032 ps 2 . 

spontaneous emission within the Markov approximation 
as being governed by a constant rate 7i 20 i 21 i 22 The mas- 
ter equation thus has the following form: 

P = -i{H cS p- pH\ s ) _ tC ^)^*) - £i(t)<?xP(rz 

- €i(t)a z p<r x - £, 2 (t)<Typa z - &(t)a z pa y ) + j<T-p<r + , 

(2) 

where H cS = [(0/2) + (&(t)/4)] a x - (£i(t)/4)<r, - 
(i/(t)/2)a z — (z7/2)<t+(T_ is non-Hermitian. 
Defining the bath correlation function C(t) — 
J °° duiJ{u)) (cos (cut) coth (f3u)/2) — i sin (tut)), the 
following time-dependent rates are obtained: 

i r* 

D(t) = — / dt'Re{C{t')}{v 2 + n 2 cos (Et')), (3a) 
& Jo 

^(*) = §j[ dt'C(t') sin (Et'), (3c) 

where E 2 = v 2 + il 2 . The environment-shifted detun- 
ing entering H c g is v (t) = v + f*dt'Im{C(t')}, and 
f3 = 1/kgT, where kg is Boltzmann's constant and T 
the temperature. The effect of the phonon bath is there- 
fore to induce frequency shifts, in u'(t) and through the 
real parts of £1(2) (t) (arising from vacuum and thermal 
fluctuations), as well as irreversible terms through the de- 
phasing D(t) and the imaginary parts of £1(2) (t) (damp- 
ing terms not present for = 0)^ 

As written above, Eq. @ is non-Markovian in the 
phonon terms, but local in time. Of primary interest 
here is the long-time behaviour of the second-order cor- 
relation function, measured under continuous excitation 
conditions. In this case, transients in the dynamics 
are ignored and the upper limits of integration in 
Eqs. (|3"a|) - (|3"c|) taken to infinity (a Markov approxi- 
mation 21 ). The rates then become time- independent, 
the key point being that they are still functions of 
the laser detuning v and dot-laser coupling il. In 
fact, only D(oo) = D and Re(^2(oo)) = Re(^) 
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will enter the final expression for g^(r), here 
given by D = (nn 2 /2E 2 ) J(E) coth ((3E/2) = 
(ann 2 E/2)e-( E /^ 2 coth (/3J5/2), and for ele- 
vated temperatures (fcsT > u c ), Re(^2) = 

&°p.[2yfFEe-(£)* {12 + £ 2 /3 2 )ern(£/u; c ) - w c (24 + 
(2E 2 + oj 2 )f3 2 )], where erfi(z) is the imaginary error func- 
tion. In addition, an expression for the bath-shifted de- 
tuning is needed; v 1 = v— J °° duJ(u))/ui = v — ay/iruj^/A. 

Both D and Re(£ 2 ) are plotted against SI in Fig. [1] for 
the resonance condition v' = 0, i.e. v = ay/nu> 3 /4. In 
particular, the dephasing rate D grows quickly as the 
driving strength is increased from zero (D oc SI 2 at small 
17) until the cut-off begins to dominate. Consequently, 
any excitonic oscillations induced by the coherent exci- 
tation become strongly damped as the driving frequency 
approaches the cut-off, as anticipated by the earlier res- 
onance arguments. For example, when 17 = lu c ~ 1 meV, 
we have 1/17 = 1.5 ps _1 , compared to D = 0.6 ps _1 . 
This damping weakens however, and eventually becomes 
negligible, as 17 increases further, beyond u> c . 

The resulting dynamics is clearly illustrated from the 
matrix form of Eq. @, (s) = M(s) — k, in terms of the 
Bloch vector s = (er_, cj + , cr z ) T , with 

- (.D + 7/2) ift'-Re(£i) 
M- I -(D+7/2) -ffl'-Re(^i) 

iQ,/2 -in/2 -7 

and k = (ilm(5 ) + Im(^), Im(^) - ilm(d), ~7) T . Here, 
17' = 17 + Re(^2) and v' = 0. The dynamic evolution of 
(°z(0) ma y now be found simply by solving for (s), while 
setting (s) = gives the required steady-state solution for 
((T + <7_). From Eq. (|TJ) we then find an analytical expres- 
sion for the normalised second-order correlation function 



gW(r) = 1 - e -<- 2D+3 ^/ 4 ('cos Ar + 2D ^ sinAr 



describing the detection statistics of photons emitted 
from the resonantly driven, dephasing QD, with A = 
■s/hlQ + Re(6)) - (2D - 7) 2 /16 now the effective oscil- 
lation frequency. The correlation function is always zero 
at t = 0, displaying the expected antibunching dip of a 
nonclassical emitter.— As r increases beyond the char- 
acteristic dephasing time (i.e. beyond min(l/D, I/7)), 
g( 2 )(r) — » 1 indicating independent photon detection due 
to dephasing of the optical transition. At intermediate 
times, the correlation function describes an exponentially 
damped oscillation provided that A is real. For a typical 
QD this is usually the case beyond weak driving fields, 
here satisfied when 17 > 1 /ieV for the parameters of 
Fig. [TJ Most importantly, damping of the induced oscil- 
lations is governed not only by 7 but also by the rate 
D, which is itself a function of the driving strength SI, 
leading to driving-dependent variations in the quality of 
oscillations as il is varied. 

This behaviour is illustrated in Fig. [5] (a) where the 
correlation function is plotted as a function of delay time 
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FIG. 2: (Color online) Top: against scaled delay time t/A 
for Q, = 2 fj,eV (red, dotted line), SI = 75 fj,eV (yellow, solid 
line), fi = 0.4 meV (blue, dashed line), and Q — 2.5 meV 
(purple, dash-dot line). Bottom: Height of the first (blue, 
solid line), second (yellow, dashed line), and third (red, dot- 
ted line) peaks above as a function of Q,. Inset: Height of 
corresponding peaks for a constant pure dephasing rate of 
(1/150) ps _1 . 7 = 1/250 ps _1 , all other parameters as in 
Fig. ID 



scaled by 1/A, to allow for easier comparison. We see 
that at the smallest driving frequency (fi = 2 fieV) no os- 
cillations are yet visible and the function exhibits only the 
antibunching dip. As the driving strength is increased 
oscillations in g^ 2 '{r) appear as a result of the induced 
coherent exciton evolution within the QD, shown here 
for SI = 75 /ieV. These oscillations are damped due 
to dephasing of the excitonic state by phonons and as 
a result of spontaneous emission, though at least three 
periods are clearly observable. However, as the driving 
strength is further increased the quality of oscillations 
declines, as seen for S7 = 0.4 meV. This is in marked 
contrast to the expectation for a constant phenomeno- 
logical pure dephasing rate r^ i 20 ' 24 where the quality 
should simply increase with SI due to a gain in Rabi fre- 
quency relative to the dephasing rate. As stated above, 
the origin of this effect is in the driving dependence of 
the dephasing rate D, itself a result of the form of the 
exciton-phonon interaction and its frequency dependence 
as seen in J(uS). For SI = 0.4 meV, D has grown suffi- 
ciently in comparison to that the oscillations are now 
strongly damped. Only when SI is increased to higher 
values, beyond the cut-off, does the damping weaken; in 
the case of SI = 2.5 meV, D is now small enough that the 
dephasing rate is dominated by 7. This point is empha- 
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FIG. 3: (Color online) Top: <?' ' against delay time r for f2 = 
75 fieV (yellow, solid line) and fi = 0.4 meV (blue, dashed 
line). Inset: As main but for constant pure dephasing rate of 
(1/150) ps -1 . Bottom: Oscillations in as a function of 
at fixed delay r = 20 ps, akin to the Rabi rotations behaviour 
in Ref. fiol All other parameters are as in Fig. [1] 

sised in Fig. [2] (b) where the heights of the three peaks 
in Fig. [2] (a) are plotted as a function of tt. After an ini- 
tial increase a dip is observed in the height of all peaks, 
and almost complete suppression in the case of the third. 
The width of this dip corresponds approximately to the 
width of the peak in D (see Fig. [lj. Of course, as f2 in- 
creases further, the peak heights rise once more reaching 
their maximum value around f2 = 3 meV. Again, this is 
contrary to the behaviour expected for a constant pure 
dephasing rate, shown in the inset for = (1/150) ps -1 , 
which exhibits no dip whatsoever. 



We can see then that for a significant range of the 
parameters considered, oscillations in (r) completely 
disappear after the first few picoseconds, resulting in in- 
dependent photodetection events after this time delay. 
This effect is shown in Fig. [3] (a) where, for Q, = 0.4 meV, 
independent detection events are expected at delays be- 
yond r ~ 40 ps, while oscillations persist over a much 
longer timescale for weaker driving of 17 = 75 fj,eV. As 
before, this can be contrasted to the case of a constant 
dephasing rate, illustrated in the inset. For both field 
strengths oscillations persist well beyond 100 ps, decay- 
ing at the same rate. Finally, in Fig [3] (b) we see that 
by fixing the time delay, here at 20 ps though the actual 
value is not crucial, the appearance, disappearance, and 
subsequent reappearance of oscillations in the correlation 
function can be observed as f2 is increased. This cor- 
responds to a crossover from correlated to uncorrelated 
photodetection and back again; a striking effect with its 
origin solely in the driving dependence of the dephasing 
rate. Besides being of fundamental interest, this signa- 
ture provides a means to map out the effective dephasing 
rate of the QD system, thus allowing for the determi- 
nation of optimum performance conditions for the dot 
operating either as a qubit or single photon source. 

To summarise, direct measurements of the second- 
order fluorescence intensity correlation function of a co- 
herently driven QD can be extremely sensitive to the 
internal dot coherence properties, going beyond what 
would be expected from more simplified models. I have 
shown that natural variations in the phonon-induced de- 
phasing rate result in regimes of correlated and uncorre- 
lated photon emission, dependent upon the driving field 
strength. This response should be observable in QDs 
by exploiting recently developed semiconductor and laser 
technology, and could also be relevant to other solid- 
state systems.— 
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